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1.  Introduction 

In  this  presentation,  we  outline  a  framework  for  Automatic  Target  Recognition  (ATR) 
in  Synthetic  Aperture  Radar  (SAR)  that  is  based  on  an  adaptive  processing  of  the 
digitally-spotlighted  phase  history  data  of  individual  targets.  In  the  conventional  SAR- 
ATR  algorithms,  the  magnitude  only  or  the  complex  SAR  image  of  a  test  target  chip  is 
analyzed  and  compared  with  a  set  of  reference  target  chips  to  determine  the  test  target 
type  or  class.  Depending  on  the  relative  coordinates  of  a  test  target  in  the  imaging  scene 
and  the  flight  path  of  the  radar-carrying  aircraft,  the  SAR  image  (chip)  of  the  test  target 
possesses  a  spatial  warping  with  respect  to  the  reference  target  chips  that  should  be 
blindly  incorporated  and/or  compensated  in  the  ATR  algorithm.  Furthermore,  the  formed 
SAR  images  exhibit  certain  (slant  plane)  parametric  variations  and  erroneous  shifts  that 
depend  on  the  type  of  the  SAR  imaging  algorithm  that  is  used.  For  the  success  of  a  SAR- 
ATR  algorithm,  the  spatial  warping  and  slant  plane  parameters  of  the  SAR  imaging 
algorithm  should  be  available  to  the  user.  In  practice,  these  parameters  of  the  formed 
SAR  image  are  not  exactly  known  to  the  SAR- ATR  user. 

Moreover,  any  comparison  between  the  test  and  reference  target  chips  requires 
understanding  and  incorporating  the  sensor  and  platform  variations  in  the  corresponding 
SAR  data  acquisitions.  The  sensor  variations  are  caused  by  various  subtle  changes 
(imperfections)  in  the  radar  system  circuitry  (e.g.,  waveform  generator,  cables,  etc.),  and 
undesirable  amplitude/phase  fluctuations  in  the  radiation  pattern  of  the  physical  radar 
between  the  reference  and  test  data  collections;  these  are  unknown  and  result  in  different 
2D  Image  Point  Response  (IPR)  or  Point  Spread  Function  (PSF)  in  the  reference  and  test 
SAR  images.  The  above-mentioned  ambiguities  result  in  unknown  subtle  geometric 
distortions  and  complex  PSF  variations  in  the  reconstructed  SAR  image  that  have  adverse 
effects  on  the  performance  of  any  SAR-ATR  algorithm. 
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In  this  work,  we  provide  a  two-stage  SAR-ATR  algorithm  that  is  not  sensitive  to  the 
above-mentioned  problems.  In  the  proposed  approach,  the  SAR-ATR  user  applies  a 
relatively  fast  CFAR  detection  algorithm  with  a  high  probability  of  false  alarm  as  well  as 
a  high  probability  of  detection  to  identify  the  coordinates  of  suspected  targets  in  a  SAR 
imaging  scene.  Next,  the  coordinates  of  each  suspected  target  is  passed  to  what  we  refer 
to  as  the  digital-spotlighting  algorithm  that  extracts  the  SAR  phase  history  data  of  a 
neighborhood  within  the  input  coordinates.  The  resultant  database  is  the  complex  SAR 
signature  (phase  history)  of  the  chip  area  at  the  desired  coordinates  that  carries 
information  on  the  variations  of  the  target's  complex  radar  signature  with  respect  to  the 
fast-time  (radar)  frequency  and  aspect  angle.  This  database  is  then  compared  with  a  set  of 
reference  phase  history  models  to  identify  the  target  at  the  input  coordinates  where  the 
original  ATR  algorithm  detected  a  suspected  target. 

A  2D  adaptive  method  for  the  above-mentioned  matching  of  test  and  reference  phase 
history  data  is  utilized,  an  algorithm  that  we  call  signal  subspace  matched  filtering,  that  is 
not  sensitive  to  the  calibration  errors  of  the  SAR  system  (i.e.,  variations  of  the  radar 
radiation  pattern  from  one  experiment  to  another).  Such  calibration  errors,  that  could 
drastically  alter  the  phase  information  of  a  target's  SAR  signature,  are  one  of  the  major 
obstacles  in  exploiting  the  complex  SAR  signature  of  a  target  in  the 
classification/recognition  problems.  The  signal  subspace  processing  method  performs  a 
blind  calibration  of  variations  of  the  IPRs  of  the  reference  and  test  SAR  phase  history 
data  using  2D  adaptive  filtering  methods.  This  calibration  also  compensates  for  variations 
in  the  IPRs  that  are  due  to  imperfect  (errors  in)  motion  data  and  required  2D  auto- 
focusing  to  compensate  for  them  (among  other  SAR  system  phase  errors  such  as  range- 
gate  slip)  that  never  results  in  the  theoretical  (ideal)  SAR  image. 


2.  Sources  of  Calibration  Errors 

SAR  signature  of  a  target  depends  on  the  radar  sensor  characteristics  (beam  pattern, 
signal  generator,  etc.)  as  well  as  the  physical  properties  of  the  target.  Figure  1  shows  the 
overall  system  diagram  for  generation  of  a  target  signature  in  SAR.  In  SAR-ATR, 
comparison  between  the  test  and  reference  target  chips  requires  understanding  and 
incorporating  the  sensor  and  platform  variations  in  the  corresponding  SAR  data 
acquisitions.  The  sensor  variations  are  caused  by  various  subtle  changes  (imperfections) 
in  the  radar  system  circuitry  (e.g.,  waveform  generator,  cables,  etc.),  and  undesirable 
amplitude/phase  fluctuations  in  the  radiation  pattern  of  the  physical  radar  between  the 
reference  and  test  data  collections;  see  Figures  2a  and  2c.  These  are  unknown  and  result 
in  different  2D  Image  Point  Response  (IPR)  or  Point  Spread  Function  (PSF)  in  the 
reference  and  test  SAR  images.  These  ambiguities  result  in  unknown  subtle  geometric 
distortions  and  complex  PSF  variations  in  the  reconstructed  SAR  image  that  have  adverse 
effects  on  the  performance  of  any  SAR-ATR  algorithm.  Variations  in  the  target  may  also 
exist;  e.g.,  rotation  in  the  gun  of  a  tank,  opening  or  closing  of  driver/commander/turret 
hatch,  etc.;  see  Figure  2d. 
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(AM  Signal) 


(PM  Signal) 


Figure  1.  Model  for  Generation  of  SAR  Signal 
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Variations  in  Radar 
Transmit-Mode 
Beam  Pattern 


Variations  in 
Motion  & 
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Transmit  Mode  Radar 
Amplitude  Pattern 
(AM  Signal) 

*T  (®>  *n*Vn-u  ) 
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Phase  Delay 
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Figure  2a.  SAR  Calibration  Errors:  Transmit  Mode 


Figure  2b.  SAR  Calibration  Errors:  Variations  in  Target 
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Amplitude  Pattern 
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Figure  2c.  SAR  Calibration  Errors:  Receive  Mode 

3.  Adaptive  Calibration  of  Dual  SAR  Imagery  via  Signal 

Subspace  Processing 

In  SAR-ATR,  the  task  of  a  signal  processor  is  to  blindly  compensate  for  the  above- 
mentioned  calibration  errors.  For  this  purpose,  we  utilize  a  2D  adaptive  method  for  the 
matching  of  test  and  reference  target  chips  is  utilized  [1].  This  algorithm,  that  we  call 
signal  subspace  matched  filtering,  is  not  sensitive  to  the  calibration  errors  of  the  SAR 
system  (i.e.,  variations  of  the  radar  radiation  pattern,  etc.  from  one  experiment  to  another) 
as  well  as  small  variations  in  the  target.  The  mathematical  foundation  of  this  approach  is 
the  same  as  an  adaptive  filtering  algorithm  that  we  have  developed  for  coherent  change 
detection  [2],  [1,  eh.  8],  and  moving  target  detection  [1,  eh.  8],  [3],  [4]  in  SAR  systems. 

The  following  outlines  the  foundation  of  the  adaptive  signal  subspace  algorithm  using 
first  a  one-dimensional  signal  space.  The  results  are  then  extended  to  the  two- 
dimensional  problems  of  SAR. 
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Signal  Subspace  Processing: 
AID  Example 


We  want  to  determine  h(x)  such  that  ffx)  =f2(x)  (i.e.  e(x)  =  0)  ,for  the  example  this  means 
that  we  must  determine  h(i)  for  i=  -2,  The  equations  for  the  output  ofh(x)  are: 


f(-2)h(-2)+  /,(-l)*(-l)+  /,(0)A(0)+  f{\m)+  f,(2)h(2)=  0) 

ff-\)h(-2)+  f(m0)+  f(2)h(l)+  f(3)h(2)=  fM(  1) 

f(0)h(-2)+  /,(l)/t(-l)  +  f(2)h(0)+  /,(3)/t(l)  +  f(4)h(2)  =  /„,(  2) 

/,  (n  ~  2)h(-2)  +/,(»-  l)/*(-l)  +  /,  (n)A(O)  +  /,(»  + 1)6(1)  +  /  (»  +  2)6(2)  =  /„,  (») 
This  can  be  represented  in  vector  form  as  shown  on  the  next  slide .... 


'  /i(-2)  " 

'  /,(-!)  ' 

7i(0)" 

'  7(1)  ' 

'  7(2)  ' 

"7,(0)' 

/(-l) 

h(~  2)  + 

/,(0) 

/»(-!)  + 

/(l) 

h{  0)  + 

/t(2) 

6(1)  + 

m 

6(2)  = 

7,0) 

1 

1 

R 

_/i  («  +  !)_ 

/,(«  +  2). 

.7(7 

V-2 


V-1 


V0  V! 


V2 


V 


out 


The  vector  representation  shows  that  the  output  resides  in  the  signal  subspace  spanned  by  V  -  {v(. ,  /'  -  -2, ...,  2}.  We  | 
can  create  an  orthonormal  basis  for  this  subspace  (via  Gram-Schmidt)  and  obtain  modified  filter  coefficients  based  on 

these  new  basis  vectors.  That  is:  I 

v_2h(- 2)  +  1)  +  v0/z(0)  +  vxh(  1)  +  v2/z(2)  = 

x_2h(-2)  +  v_jA(-l)  +  v0/i(0)  +  v^l)  +  v2/z(2) 


Where  the  notation  y .  indicates  that  the  vector  is  part  of  the  orthonormal  basis  and  the  notation  h(i)  indicates  that  | 
the  coefficient  has  been  adjusted  appropriately  (for  the  orthonormal  basis  vector). 
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We  can  solve  for  each  of  the  filter  coefficients  by  calculating  the  inner  product  with  the  corresponding  orthonormal 

basis  vector: 


h(i)  =  (x_2h{-2)  +  v_j/z(— 1)  +  \0h(0)  +  VjA(l)  +  \2h(  2), 

=  (v/,vMI/) 


These  coefficients  can  then  be  used  to  calculate  the  projection  of f2(x)  signal  into  the  signal  subspace  of f/x). 


^  The  resulting  output  of  h(x)  (the  estimate  of f2(x))  is  simply  the  projection  of f2(x)  into  the 

signal  subspace  of f/x)  !! 


The  difference  between  the  reference  signal  and  the  projected  test  signal  becomes:  e  =  f2—  fx 
which  also  shows  how  much  the  test  signal  differs  from  the  reference. 
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Extension  to  2D 


e(x,y) 


Linear  Filter: 

f2k)(Xi’yj)=  Z  S  hmn  flk)(Xi-mAx’yj-nAy) 

m=—nx  n--ny 

Signal  Subspace: 

=  [f(k\x.  - mAx,y.  - nAy);  m  =  -nx,---,nx  ,n=-ny,---,ny] 


The  orthonormalized  Signal  Subspace: 

®[k)  =  \d^(xl -mAx,y  y  -nAy);  m  =  -nx,---,nx  ,  n  =  -ny,---,ny 

X.  =  X  .  r  ,  X  .  r  r  , 

i  mm  ,k  mm  ,k  max,  a: 


y  .  =  y  •  , ,  y  •  1  +1  ,...,y  1  ] 

j  min  ,k  mm,  a:  max,A:J 


The  filter  output  (for  the  normalized  subspace): 

flk)  (Xi  ,y j)=  E  Z  flk\Xi  ~  mA<  >  -  « A>. ) 


m=-nx  n=-nv 


Analogous  to  the  ID  case,  the  output  of  the  filter  is  simply  the 
projection  of f2<k>(x,y)  into  @(k). 
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4.  Adaptive  SAR-ATR:  Ka  Band  Turntable 

(ISAR)  Data 


We  next  examine  the  application  of  the  SSP  algorithm  in  SAR-ATR  using  a  Ka  band 
turntable  ISAR  database.  Figure  3a  shows  a  portion  of  the  measured  ISAR  data  in  the 
radar  frequency  and  aspect  angle  domain.  The  first  issue  that  we  notice  for  this  database 
is  the  irregular  vertical  lines.  Figure  3b  shows  the  cumulative  spectrum  of  the  data.  This 
indicates  irregular  transmitter  power  variations  in  the  ISAR  data. 

This  is  also  a  calibration  error  source  that  can  be  easily  compensated  for  using  the 
inverse  of  the  distribution  in  Figure  3b  as  a  filter  to  be  applied  to  the  measured  ISAR 
data.  (Note  that  this  is  magnitude  only  calibration.  One  may  also  perform  phase 
calibration  provided  an  omni-directional  target  had  been  put  in  the  imaging  scene.) 
Figure  3c  shows  the  magnitude-calibrated  ISAR  data. 


Measured  SAR  Data 


Radar  Frequency,  GFIz 


Figure  3a.  Measured  ISAR  Data 
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Cumulative  Relative  Magnitude 


Figure  3b.  Cumulative  Spectrum  of  Measured  ISAR  Data 


We  use  two  SAR/ISAR  reconstruction  algorithms  to  form  images  from  the  magnitude- 
calibrated  data.  The  first  method  is  the  polar  format  algorithm  that  is  based  on 
approximations;  the  second  method  is  the  wavefront  reconstruction  algorithm  that  is 
error-free  [1,  ch.  7].  Figures  4a  and  4b  are  the  resultant  images  using  16  degrees  of 
integration  angle;  similar  results  for  a  32-degree  integration  angle  are  shown  in  Figures 
5a  and  5b.  Note  that  the  wavefront  reconstruction  algorithm  provides  better  images.  For 
our  SAR- ATR  study,  the  waevfront  images  are  used. 

For  the  SAR-ATR  study,  we  consider  two  of  the  ISAR  databases  of  a  T-72  tank: 
T06.frq  and  T07.frq.  The  difference  between  the  two  databases  is  the  orientation  of  the 
gun  of  the  vehicle:  for  the  T06.frq  data,  the  gun  is  at  the  azimuth  angle  of  -20  degrees;  the 
azimuth  angle  of  the  gun  is  zero  for  the  T07.ftq  data.  Thus,  the  recognition  algorithm 
should  match  the  structure  of  the  T-72  tank,  except  for  the  gun,  in  the  two  databases. 


21  -10 


RTO-MP-SET-096 


UNCLASSIFIED/UNLIMITED 


NATO 

_ 1 _ 

OTAN 

UNCLASSIFIED/UNLIMITED 

Adaptive  Processing  of  SAR  Data  for  ATR 


Calibarted  SAR  Data 


33.7  33.8  33.9  34  34.1  34.2  34.3 

Radar  Frequency,  GFIz 


Figure  3c.  Calibrated  ISAR  Data 


A  simple  differencing  or  correlation  of  the  two  ISAR  images  would  not  indicate  a 
match.  Figures  6a  and  6b,  respectively,  show  the  incoming  and  outgoing  signal  subspace 
difference  (SSD)  images.  Note  that  the  T-72  components  are  matched  except  for  the  gun 
and  an  area  near  the  upper  hatch.  The  latter  might  be  an  unintended  change  that  occurred 
while  moving  the  tank  on  and  off  the  turntable.  This  can  also  be  seen  in  Figure  7a  and  7b 
that  are  the  close  ups  of  the  original  ISAR  images  and  SSD  images. 
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Figure  4a.  Polar  Format  Reconstruction: 
Integration  Angle  =  16  Degrees 
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Figure  4b.  Wavefront  Reconstruction: 
Integration  Angle  =  16  Degrees 


RTO-MP-SET-096 


21  -13 


UNCLASSIFIED/UNLIMITED 


UNCLASSIFIED/UNLIMITED 


Adaptive  Processing  of  SAR  Data  for  ATR 


ORGANIZATION 


Figure  5a.  Polar  Format  Reconstruction: 
Integration  Angle  =  32  Degrees 
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Figure  5b.  Wavefront  Reconstruction: 
Integration  Angle  =  32  Degrees 
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Incoming  Coherent  SSD 
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Figure  6a.  Incoming  Coherent  Signal  Subspace 

Difference 
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Outgoing  Coherent  SSD 
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Figure  6b.  Outgoing  Coherent  Signal  Subspace 

Difference 
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Reference  Image:  TOG 
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Figure  7a 
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ORGANIZATION 


5.  Summary 


Sensor  calibration  is  a  critical  problem  that  has  to  be  dealt  with  (though  ignored)  in 
noncoherent  and  coherent  imaging  systems  such  as  SAR.  Simple  modification  and/or 
reformulation  of  existing  ID  adaptive  (blind)  filtering  methods  exist  for  2D  calibration 
problems  of  SAR-ATR,  SAR-Coherent  Change  Detection  (CCD),  SAR-MTI,  etc.  This 
paper  exhibited  an  example  of  the  application  this  method  in  the  SAR-ATR  problem. 
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Sources  of  Calibration 
Errors  in  SAR 


Sources  of  Calibration  Errors 


•  SAR  signature  of  a  target  depends  on  the  radar 
sensor  characteristics  (beam  pattern,  signal 
generator,  etc.)  as  well  as  the  physical  properties 
of  the  target 

•  In  SAR-ATR,  comparison  between  the  test  and 
reference  target  chips  requires  understanding 
and  incorporating  the  sensor  and  platform 
variations  in  the  corresponding  SAR  data 
acquisitions 


Sources  of  Calibration  Errors, 

Cont. 


The  sensor  variations  are  caused  by  various  subtle 
changes  (imperfections)  in  the  radar  system  circuitry 
(e.g.,  waveform  generator,  cables,  etc.),  and 
undesirable  amplitude/phase  fluctuations  in  the 
radiation  pattern  of  the  physical  radar  between  the 
reference  and  test  data  collections 

These  are  unknown  and  result  in  different  2D  Image 
Point  Response  (IPR)  or  Point  Spread  Function  (PSF) 
in  the  reference  and  test  SAR  images 


Sources  of  Calibration  Errors, 

Cont. 


These  ambiguities  result  in  unknown  subtle 
geometric  distortions  and  complex  PSF  variations  in 
the  reconstructed  SAR  image  that  have  adverse 
effects  on  the  performance  of  any  SAR-ATR 
algorithm 

Variations  in  the  target  may  also  exist;  e.g.,  rotation 
in  the  gun  of  a  tank,  opening  or  closing  of 
driver/commander/turret  hatch,  etc. 


Model  for  Generation 
of  SAR  Signal 


Transmit  Mode  Radar 
Amplitude  Pattern 
(AM  Signal) 

aT  <w,*„.y„-u} 


Overall 
SAR  System 
Signal  Model 

r  " 


Receive  Mode  Radar 
Amplitude  Pattern 
(AM  Signal) 


Radar-tO-Target 
Phase  Delay 
(PM  Signal) 

e*p  [  -jk  \j  7}  +  (  yn-  u  ) 2  ] 


Target 

Amplitude  Pattern 
(AM  Signal) 

a„  I".  xn.  yr-u  ) 


Target-to- Radar 
Phase  Delay 
(PM  Signal) 


Variations  in  Radar 
Transmit-Mode 
Beam  Pattern 


Transmit  Mods  Radar 
Amplitude  Pattern 
(AM  Signal) 

*T  (w.  *n*yn'u  ) 


exp 


1 


Variations  in 
Motion  & 
Autofocus 

Radar-to-Target 
Phase  Delay 
(PM  Signal) 

[-ik  \j  *„2  *(y„-u)2 

X) - 1 


Variations  in  Target 
Coherent  SAR 
Signature 


Radar-Target 
Radiation  Pattern 


Target 

Amplitude  Pattern 
(AM  Signal} 

an  (<*>  *„>  V,Tl<  ) 


(AM  Signal)  (PM  Signal) 


Variations  in  Radar 
Receive-Mode 
Beam  Pattern 


Variations  in 
Motion  & 
Autofocus 


Adaptive  Calibration  of 
Dual  SAR  Imagery  via 
Signal  Subspace 
Processing 


Adaptive  SAR-ATR  Solution 


A  2D  adaptive  method  for  the  matching  of  test  and 
reference  target  chips  is  utilized 

This  algorithm,  that  we  call  signal  subspace  matched 
filtering,  is  not  sensitive  to  the  calibration  errors  of  the 
SAR  system  (i.e.,  variations  of  the  radar  radiation 
pattern,  etc.  from  one  experiment  to  another)  as  well 
as  small  variations  in  the  target 


Signal  Subspace  Processing: 

AID  Example 


We  want  to  determine  h(x)  such  that  f(x)  —f2(x)  (i.e.  e(x)  =  0)  ,for  the  example  this  means 
that  we  must  determine  h(i)  for  i=  The  equations  for  the  output  of  h(x)  are: 


f(-2)h(-2)+  f(- 1)*(-1)+  /(0)*(0)+  /,(l)/t(l)+  f(2)h(2)=  0) 

f{-\)h{-2)+  f(0)h(-l)+  f(l)h(0)+  f(2)h(l)+  f(3)h(2)=  /„„,(!) 

Ami-2)  +  f(l)h(-l)+  f(2)h(0)+ f(3)h(l)+  f(4)h(2)=  fM(  2) 


/i (n ~  2)h(-2)  +  fx{n- +  fx (n)h( 0)  +  fx(n  +  \)h{\)  +  fx(n  +  2)h(2)  =  fout ( n ) 


This  can  be  represented  in  vector  form  as  shown  on  the  next  slide. . .. 


Signal  Subspace  Processing 
AID  Example,  Cont. 
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The  vector  representation  shows  that  the  output  resides  in  the  signal  subspace  spanned  by  V  =  {v. ,  i  =  -2,  . . 2} .  We 
can  create  an  orthonormal  basis  for  this  subspace  (via  Gram-Schmidt)  and  obtain  modified  filter  coefficients  based  on 

these  new  basis  vectors.  That  is: 


\_2h(-  2)  +  1)  +  v0/*(0)  +  v^l)  +  v2/i(2)  = 

y_2h{-2)  +  v_|l( — 1)  +  \0h(0)  +  Vj^(l)  +  v2/i(2) 


Where  the  notation  y  indicates  that  the  vector  is  part  of  the  orthonormal  basis  and  the  notation  h(i)  indicates  that 
the  coefficient  has  been  adjusted  appropriately  (for  the  orthonormal  basis  vector). 


Signal  Subspace  Processing: 
AID  Example,  Cont. 


We  can  solve  for  each  of  the  filter  coefficients  by  calculating  the  inner  product  with  the  corresponding  orthonormal 

basis  vector: 


h(i)  -  (\_2h(- 2)  +  v_j^(-l)  +  v0/?(0)  +  v1/z(l)  +  v2h( 2),  v; 

=  (Vi,Vout) 


These  coefficients  can  then  be  used  to  calculate  the  projection  of f2(x)  signal  into  the  signal  subspace  of fj(x). 


4  The  resulting  output  of  h(x)  (the  estimate  of f2(x))  is  simply  the  projection  of f2(x)  into  the 

signal  subspace  of 'f2(x)  ! ! 


/V 

The  difference  between  the  reference  signal  and  the  projected  test  signal  becomes:  e  =  f2~  f\ 
which  also  shows  how  much  the  test  signal  differs  from  the  reference. 


Extension  to  2D 


Linear  Filter: 

flk)  (*/  >yj)=  Z  Z  hmn  flk)  (Xi  ~  mAx  >  y}  -  nAy  ) 

m=-nx  n=-ny 

Signal  Subspace: 

(jjW  =[/«(*. ,-mAx,y. -nAy);  m  =  -nx,  —  ,nx  ,n  =  -ny,---,ny] 


Extension  to  2D,  Cont 


The  orthonormalized  Signal  Subspace: 

©!*’  =W\x,-mkx,y  j-nA  );  m  =  -nx,--,nx  ,  n  = 


X.  =  X  .  7  ,  X  .  7  +1,...  ,x 

i  min,  A:  min,  A:  max.  A: 

y  .  =  y  •  7 ,  y  •  7  +1  ,...,y  7  ] 

'y  ymin,A;  ymin,A:  vma  x,A:J 


The  filter  output  (for  the  normalized  subspace): 

fik} (X, ,y})=  X  X  h(y  yk) (xt  -mAx,yj  -nAy) 


m--nv  n=—n. 


Analogous  to  the  ID  case,  the  output  of  the  filter  is  simply  the 
projection  of f2<k>(x,y)  into  ©(k). 


Signal  Subspace  Matching  or 
Change  Detection 
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Adaptive  SAR-ATR: 
Ka  Band  Turntable 
(ISAR)  Data 
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Wavefront  CSAR  Reconstruction:  Integration  Angle  =  16 
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Wavefront  CSAR  Reconstruction:  Integration  Angle  =  32 
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Block  Partitioning  of  Reconstruction 
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Test  Image:  T07 
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Summary 


Sensor  calibration  is  a  critical  problem  that 
has  to  be  dealt  with  (though  ignored)  in 
noncoherent  and  coherent  imaging 
systems  such  as  SAR 

Simple  modification  and/or  reformulation 
of  existing  ID  adaptive  (blind)  filtering 
methods  exist  for  2D  calibration  problems 
of  SAR-ATR,  SAR-Coherent  Change 
Detection  (CCD),  SAR-MTI,  etc. 
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